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A  kinetic  model  of  a  porous  Pt-Ru  anode  for  methanol  oxidation  is  presented.  It  was  based  on  the  dual¬ 
site  mechanism  for  methanol  oxidation  and  used  to  predict  anode  performance  and  the  influence  of 
species  adsorption  on  the  overall  oxidation  (macro-)  kinetics.  The  performance  of  the  porous  Pt-Ru 
anode  depended  on  the  parameters  of  the  intrinsic  chemical  kinetics  of  methanol  oxidation  and  physical 
parameters  such  as  electrode  thickness,  surface  area,  effective  diffusion  and  charge  transfer  coefficients 
and  concentration  of  methanol  and  temperature.  The  model  was  solved  by  using  the  finite  difference 
method  with  a  subroutine  for  solving  a  set  of  nonlinear  algebraic  equations  in  each  step.  Surface  cover¬ 
age  ratio  distributions  of  adsorbed  species,  effectiveness  of  the  porous  electrode  and  macro-polarisation 
curves  were  obtained.  The  simulated  polarisation  curves  were  compared  to  experimental  polarisation 
data  for  methanol  oxidation  on  Pt-Ru  porous  anodes  at  different  temperatures  and  methanol  concen¬ 
trations.  The  intrinsic  kinetic  parameters  were  regressed  from  the  corresponding  experimental  data.  The 
predicted  polarisation  curves  calculated  by  the  model,  were  consistent  with  experimental  polarisation 
data  at  lower  current  densities.  The  departure  of  experimental  data  from  the  predicted  polarisation  curves 
at  high  concentration  and  high  apparent  current  densities  was  believed  to  be  due  to  two-phase  flow  in 
the  electrode. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  direct  methanol  fuel  cell  (DMFC)  has  a  lower  power  density, 
a  larger  voltage  loss  and  lower  limiting  current  than  the  correspond¬ 
ing  hydrogen  polymer  electrolyte  membrane  fuel  cells  (PEMFCs). 
There  are  considerable  differences  in  reported  DMFC  performances ; 
since  issues  of  catalyst,  membrane,  kinetics  and  transport  processes 
arising  with  the  DMFC  are  more  complex  than  those  from  hydrogen 
PEMFCs  [1-3]. 

A  typical  PEMFC  is  made  up  of  a  membrane  electrode  assem¬ 
bly  (MEA),  consisting  of  the  membrane,  electrocatalyst  layers  and 
reactant  diffusion  layers,  and  flow  channels  incorporating  cur¬ 
rent  collectors.  Important  components  for  cell  performance  are 
the  porous  electrodes  in  the  MEA  where  electrochemical  reactions 
occur  (see  Fig.  1).  Unlike,  hydrogen  PEMFCs,  a  dominating  kinetic 
resistance  (polarisation)  in  the  DMFC  is  at  the  anode  due  to  sluggish 
oxidation  of  methanol  [1-4]. 

Generally  porous  electrode  structures  are  used  to  support  high 
current  densities  at  low  polarisation,  although  typically  current  and 
overpotential  are  non-uniform  due  to  the  influence  of  mass  and 
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charge  transfer  in  the  structure  [5-7].  Generalized  mathematical 
models  of  porous;  three-dimensional  electrodes  have  been  devel¬ 
oped  in  previous  papers  [7-11  ].  For  isothermal  systems  the  models 
describe  the  coupled  potential  and  concentration  distributions  in 
the  porous  or  packed  bed  electrodes  at  steady  state. 

Even  though  the  mechanism  of  methanol  electro-oxidation  on 
Pt-Ru  catalyst  in  the  DMFC  is  not  fully  understood,  researchers  have 
come  to  a  common  understanding,  namely  that  simple  Tafel  or  B-V 
equations  are  unable  to  represent  the  intrinsic  methanol  oxidation 
kinetics  and  a  dual-site  mechanism  on  Pt-Ru  catalyst  surfaces  has 
been  accepted  [12-16].  Consequently  for  a  model  of  the  observed 
or  macro-kinetics  of  the  DMFC  anode  it  is  important  to  include  the 
appropriate  dual-site  mechanism  of  methanol  oxidation. 

A  few  models  have  been  developed  for  unit  cells  or  the  anode 
of  the  DMFC.  Meyers  and  Newman  developed  a  comprehensive 
model  which  describes  the  thermodynamics,  transport  phenom¬ 
ena,  and  electrode  kinetics  of  the  system  in  DMFC  [4].  Garcia 
et  al.  presented  a  semi-analytical  model  which  can  be  solved 
rapidly  so  that  it  is  suitable  for  inclusion  in  real-time  system 
level  DMFC  simulations  [17].  Nordlund  and  Lindbergh  adopted  a 
porous  model  with  agglomerates  and  kinetic  equations  based  on 
surface  coverage  to  study  the  influence  of  porous  structure  on 
the  DMFC  anode  [18].  Scott  and  Argyropoulos  developed  a  one¬ 
dimensional  model  of  the  current  distribution  in  the  anode  of 
DMFC,  which  can  be  solved  analytically  [19].  However,  combina- 


0378-7753 /$  -  see  front  matter  ©  2009  Elsevier  B.V.  All  rights  reserved. 
doi:10.1016/j.jpowsour.2009.07.028 


2 


Y.-P.  Sun  et  al.  /  Journal  of  Power  Sources  195  (2010)  1-10 


Nomenclature 

a  specific  electrode  area  (m-1 ) 

C  dimensionless  concentration 

c  concentration  (mol  m-3 ) 

D0  diffusion  coefficient  (m2  s-1 ) 

De  effective  diffusion  coefficient  (m2  s-1 ) 

F  overpotential  (V) 

F  Faraday  constant  (C  mol-1 ) 

J  dimensionless  current  density 

jai  dimensionless  adsorption  limiting  current  density 

j  current  density  (A  m-2 ) 

jo  exchange  current  density  (A  m-2 ) 

jal  adsorption  limiting  current  density  (A  m-2 ) 

k 4  surface  rate  constant  (mol  m-2  s-1 ) 

l  thickness  of  porous  anode  (m) 

Nm  flux  of  methanol  (mol  m-2  s-1 ) 

n  number  of  electrons  involved  in  reaction 

Rg  universal  gas  constant  (J  mol-1 I<-1 ) 

R  apparent  reaction  rate  (mol  m-2  s-1 ) 

r4  surface  reaction  rate  (mol  m-2  s-1 ) 

s  dimensionless  modulus,  aj0l2 l{nFDec^) 

T  temperature  (I<) 

X  dimensionless  distance 

x  distance  (m) 

Greek  letters 

P  V-P3)FIRT=(\-P4)FIRT(V-') 

s  porosity 

0  coverage  ratio 

Ke  effective  conductivity  (£2-1  m-1 ) 

fi  dimensionless  modulus,  v2/<2>° 

v2  dimensionless  parameter,  aj0l2  fi  I  Ke 

§  effectiveness 

0°  potential  of  open-circuit  (V) 

cj)1  potential  of  electrolyte  (V) 

0m  potential  of  matrix  (V) 

dimensionless  relative  overpotential,  <2>//3F° 

Subscripts 
ads  adsorbed 

e  effective 

i  species  i 

loc  local  state 

M  methanol 

T  total  or  apparent  state 

Superscripts 

0  at  the  boundary  of  the  diffusion  layer  X  =  1 


tion  of  the  dual-site  methanol  oxidation  model  with  the  material 
and  charge  transport  behaviour  of  the  porous  DMFC  anode  has  not 
been  addressed. 

In  the  present  model,  the  dual-site  methanol  oxidation  mech¬ 
anism  [15,16]  was  combined  with  material  and  charge  balances 
[7-11]  to  simulate  the  porous  anode  behaviour  and  predict  the 
effect  of  methanol  concentration  and  temperature  on  anode  polar¬ 
isation;  thereby  providing  a  useful  tool  for  theoretical  analysis  of 
macro-kinetics. 


C02  +  H2O 


CH30H  +  H2O 


Porosity  0.3 
Thickness  =  3.8  -  50  jun 
0  _Spccific  area  =»I.O*  lO’rrT1 


FLOW 

CHANNEL 


DIFFUSION 

LAYER 


MEMBRANE 


CATALYST 

LAYER 


Effective  conductivity  =  3.4  £T*  m'1 
Effective  diffusion  coefficients 
=  2.63  *  10-'°  -  9.91  x  10 10  m2  s' 


Fig.  1.  Scheme  of  the  porous  anode  in  DMFC. 


2.  Model  equations  for  macro-kinetics 

Carbon  dioxide  is  produced  from  methanol  electrochemical  oxi¬ 
dation  at  the  anode  surface  in  acidic  solution,  in  the  DMFC  [1,20,21  ]. 

CH3OH  +  H20  -*  C02  +  6H+  4-  6e"  (1 ) 

The  formation  of  C02  bubble  is  related  to  factors,  such  as  local 
reaction  conditions,  solubility,  surface  tension,  pore  structure,  etc. 
and  may  raise  issues  associated  with  transport  phenomena  of  two- 
phase  flow.  By  choosing  appropriate  conditions,  it  is  possible  for 
nucleation  of  C02  to  take  place  in  the  diffusion  layer  away  from 
the  active  layer  of  porous  anode  or  to  be  limited  to  a  small  region 
of  the  active  layer  [20,21  ],  which  will  be  favourable  for  anode  per¬ 
formance.  In  this  model,  the  influence  of  C02  bubbling  is  omitted, 
although  it  is  indeed  difficult  to  avoid  the  effect  of  two-phase  flow 
in  the  diffusion  layer  and  flow  channel  to  analyse  the  whole  DMFC 
(Fig.l). 

The  main  task  in  the  mathematic  model  of  the  porous  anode  was 
to  introduce  the  intrinsic  kinetic  expression  for  methanol  oxidation 
into  the  charge  and  material  balances;  which  are  described  below. 

2.1.  Intrinsic  kinetics  expression 

The  intrinsic  kinetics  of  methanol  electro-oxidation  on  Pt-Ru 
catalyst  has  been  studied  in  detail  [15].  The  assumed  reaction  steps 
(from  (A.1)-(A.4))  and  the  corresponding  rate  expressions  (from 
(R.1)-(R.4))  are  listed  in  Table  1.  There  is  a  series  of  reactions  for 
key  components  of  methanol  (M).  At  a  steady  state,  the  rates  of 
these  reaction  steps  should  be  equal,  i.e.  r\  =r2  =  r4,  and  only  one 
is  needed  to  represent  the  total  reaction  rate  of  the  series.  In  this 
paper,  the  expression  for  r4  for  the  last  step  in  the  series,  involv¬ 
ing  information  of  all  elementary  steps  (including  the  parallel  OH- 
adsorption),  was  selected  as  the  intrinsic  kinetics  expression.  The 
last  step  is  a  surface  reaction  which  is  believed  to  occur  on  both 
the  Pt  and  Ru  sites  depending  upon  the  hydroxyl  ion.  Two  cases 
of  kinetics,  have  been  developed  previously  according  to  the  dual¬ 
site  mechanism,  depending  upon  whether  or  not  surface  coverage 
of  OH  on  Pt  site  could  be  neglected  [14-16].  Since  the  modeling 
procedures  are  similar  for  both  cases,  OHads  Ru  was  mainly  used  for 
modeling  in  this  paper. 

The  surface  rate  expression  (molm-2s-1)  of  the  last  reaction 
step  (Eq.  (A.4)  in  Table  1 )  can  be  written  as: 


-r4  =  /<40oh,ru#co  exp(/3F)  (2) 

where  /<4  is  the  rate  constant,  0Oh,ru  and  0C o  are  the  coverage 
ratios  of  OHads  at  Ru  sites  and  COads  at  Pt  sites  respectively, 
p  =  {\ -p3)FlRT={I  -£4)F/RT  [15,16],  and  F  is  the  overpotential. 
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Table  1 

The  reaction  steps  for  methanol  oxidation  and  their  rate  expressions  and  surface  coverage  of  intermediate  species  according  to  the  dual-site  mechanism  on  a  Pt-Ru  anode 
[15]. 


Elemental  steps 

Rate  expressions 

CH3OH&CH3OHads 

(A.l) 

— ri  =  ki  Cm(  1  -  #oH,pt  -  $co  -  #m)  -  -  I<20m  e(“2F£/,Rr) 

(R.1) 

CH3OHads^COods  +  4H+  +  4e- 
H20&OHads,pt  +  H++e- 

k' 

(A.2) 

(A.3.1) 

— r2  =  k2e M  _  k4A0OH,Pt0co  e^-WW  _  /<4,2£OH,Ru0c o 

-r3;i  =  /c3,iQh2o(1  -  #0H,pt  -  Oco  -  0m)  e((1~^3)F£/RT)  -  k'3  ^oH.pt^co  e((~^3)F£/RT) 

(R.2) 

(R.3.1) 

.3,1 

H20  420Hads.Ru  +  H+  +  e- 

k3  2 

(A.3.2) 

-r3,2  =  /<3,2aH2o(l  -  0oh,ru)  e^-^3)^)  -  k'3  20OH,RU0Co  e«-WW 

(R.3.2) 

COads  +  (OHads.Pt  +  OHadSiRu)%'C02  +  H+ 

+  e-(A.4) 

-r4  =  fe4,i0oH,Pt^co  +  k4,20OH,Ru0co 

(R.4) 

k4,2 

Surface  coverage  of  different  intermediate  species  in  two  different  cases 


Case  A  (#oH,pt  =  0) 


Case  B  (#oH,pt  #  0) 


A  a  (k\° 

0CO  =  1  -0M—  I  I  -  [ - - ) 

r\  k-iOiA  e(“2F£/RgTJ 

&CO  = 


Oco  = 


k20M  e(a2FE/RgC 

k3,2  e((l^3)F£/£gF)-k3i20oH,Ru  20OH,Ru  e«-P3)FW) 

- *4, 2%H.RU  eUl-/54)F£/«gJJ  ’ - 


_  /  k\  °M  \  _  /  k2e M  e^2F£/Rgr)  \ 

1  l  fclcM  J  V  klcM  / 

k20M  e(“2££/£gT) 

k4,iWt  e(('-/34)W£gJ)+k4j20OHRu  e((i-/W*gU 

e«^3)FW)lk'  i0OH,Pt  e((^3)££/£g£) 
- k3> !  ell  1  -  ^3  )££/£g  1 ) +fe4; " 0OH  pt  ’ eU  1  -/?4  Wg 1 T - 


#CO  =  1  -  #OH,Pt  -  #1 
#CO  = 

#CO  = 


#CO  = 


fe3,2(l-^OH.Ru)e«1-^3)££/£gF)-k'  20OH.Ru  e«^3 )F£/£gF) 
/<4,2%H,Ru  e«M4>W) 


Note:  cM  =  cOCm,£  =  £°‘//. 


The  surface  current  density  (Am-2)  will  be 

j  =  —  (-*4)  =  — /<4^oh,ru^co  exp(/3£)  (3) 

where  s*  is  the  stoichiometric  coefficient  in  the  electrode  reaction, 
and  n  is  the  number  of  electrons  involved  in  a  reaction. 

With  Sj  =  1,  we  have 

j  =  Moh.ru^co  exp(^£)  (4) 

where  jo  =  nFk4. 

Eq.  (4)  describes  the  behaviour  of  activation  polarisation  at  the 
interfacial  surface  of  pores  in  the  anode. 

It  should  be  noted  that  the  expressions  of  0Oh,ru  and  0C o  are 
complex,  which  meant  that  the  intrinsic  kinetics  expression  used 
to  analyse  the  Pt-Ru  porous  anode  was  nonlinear.  To  obtain  the 
values  of  0Oh,ru  and  0C o  a  set  of  nonlinear  algebraic  equations  [15] 
listed  in  Table  1  were  solved.  Two  cases  for  adsorption  of  OH  are 
presented  in  the  table. 

2.2.  Mass  and  charge  balances 

Transfer  processes  in  a  macro-reaction  system  can  form  spa¬ 
tial  distributions  of  reaction  kinetic  parameters,  namely,  the  mass 
and  energy  fields  [6,10].  The  fields  in  a  three-dimensional  electrode 
can  be  described  mathematically  by  mass  and  energy  balances  in  a 
differential  volume  of  the  system  [7-11  ]. 

V.J %  =  ar,  =  -2LqKE,Q)  (5) 

ntt 

v  -j  =  aj(E,  q)  =  ^—ar,  (6) 

SZ 

where  N,-  is  the  flux  density  of  species  i  (mol  m-2  s-1 ),  V  •  N,  is  the 
divergence  of  Nz  (molm-3  s-1),  V  j  is  the  divergence  of  j  (Am-3), 
Q  and  r,-  are  the  concentration  of  species  i  (molm-3)  and  rate  of 
reactant  of  species  i  (mol  m-2  s-1 )  respectively  and  a  is  the  specific 
area  of  the  electrode  (m-1 ). 

The  source  term  j(E,q)  in  Eqs.  (5)  and  (6)  is  determined  by  the 
reaction  mechanism.  For  the  porous  anode  of  the  DMFC,  when  there 
is  excess  water,  the  mass  balance  will  be  only  for  the  key  compo¬ 
nent  methanol  as  the  influence  of  mass  transfer  of  products  can  be 
neglected.  Also  the  energy  balance  will  only  consider  charge  and  not 
heat  because  an  isothermal  electrochemical  system  is  assumed.  A 


key  point  of  this  work  was  to  use  the  intrinsic  kinetic  expression  Eq. 
(2)  as  the  chemical  source  term  for  the  mass  and  charge  balances 
[7,10,11]. 

2.2.1.  Mass  balance  in  a  differential  volume 

We  assume  diffusion  as  the  only  mass  transport  mechanism  of 
methanol  in  the  porous  structure,  which  is  characterised  by  the 
effective  diffusion  coefficient  De  (m2  s-1 );  according  to  Fields  Law, 

Nm  =  -DeVcM  (7) 

where  vcM  is  the  concentration  gradient  of  methanol. 

For  a  one-dimensional  slab  of  porous  electrode,  we  substitute 
Eq.  (7),  Eq.  (2)  or  Eq.  (4)  into  Eq.  (5), 

^  =  ^-^4^oh,ru^co  exp(/3£)  =  ^^-^oh.ru^co  exp(^E)  (8) 

which  describes  the  effect  of  concentration  changes  on  activation 
polarisation  in  the  porous  anode. 

Assuming  the  flux  of  methanol  through  the  membrane  can  be 
neglected  and  that  both  diffusion  in  the  diffusion  layer  and  two- 
phase  flow  in  the  flow  channel  do  not  influence  the  reaction  rate, 
the  concentration  of  methanol  was  constant  at  the  boundary 
between  the  diffusion  layer  and  active  layer.  The  thickness  of  the 
porous  active  layer  was,  /,  then  the  boundary  conditions  for  the 
second  order  differential  equation  was 

x  =  0,  =  0  (at  membrane); 

dx 

x  =  l,  cM  =  cjJj  (at  the  boundary  of  the  diffusion  layer) 

Therefore  the  flux  density  of  methanol  into  the  active  layer,  the 
macro-reaction  rate  R,  was 

""-'--“•(tfL  (10) 


and  the  flux  density  of  charge,  the  total  current  density,  was 


2.2.2.  Charge  balance  in  a  differential  volume 

Assuming  the  ohmic  resistance  of  the  electrode  matrix  can  be 
neglected,  charge  transport  in  the  electrolyte  of  the  porous  struc- 
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ture,  which  was  characterised  by  the  effectiveness  conductivity  /ce 
(£2_1  m-1 ),  was  given  by  Ohm’s  law: 

}=-KeV0'  (12) 

where  4 1)1  is  the  potential  of  the  electrolyte  phase  (V),  v0*  is  the 
gradient  of  potential. 

The  overpotential,  E  is  written  as  E  =  0m-0/ -0°  [18],  where 
both.  (pm  is  the  potential  of  the  solid  electrode  phase  and,  0°  is 
the  open-circuit  value  can  be  considered  as  constant.  Substituting 
Eq.  (12)  into  Eq.  (6),  the  differential  charge  balance  to  describe  the 
potential  field  in  the  porous  anode  becomes  a  nonlinear  Poisson 
equation 

V2E  =  — j  (13) 

Ke 

For  a  slab  of  porous  electrode,  we  deploy  the  Laplacian  in  Eq. 
(13)  and  substitute  Eq.  (2)  or  (4)  into  Eq.  (13), 

=  fnFk4ffOHiRu0Co  exp(/3E)  =  ^6>Oh,ru0co  exP(£E)  04) 

CLX  /Ce  Ate 

which  describes  the  effect  of  the  ionic  resistance  on  activation 

polarisation  in  the  porous  anode. 

The  overpotential  is  assumed  constant  E°  at  the  boundary 
between  the  diffusion  layer  and  active  layer,  and  the  boundary 
conditions  of  the  second  order  differential  equation  will  be 

dE 

x  =  0,  —  =  0  (at  membrane) 

x  =  l,  E  =  E°  (at  the  border  of  diffusion  layer) 


Therefore  the  charge  flux  density  coming  from  the  active  layer, 
i.e.  the  total  current  density  jT,  is 


and  the  flux  density  of  methanol  going  into  the  active  layer,  the 
macro-reaction  rate  R,  is 


Coupling  Eqs.  (8)  and  (9)  with  Eqs.  (14)  and  (15),  we  obtain  a 
set  of  nonlinear  differential  equations  which  describes  the  cou¬ 
pled  concentration  and  overpotential  distributions  along  the  lateral 
direction  x  of  the  porous  electrode,  i.e.  the  mathematic  model  of 
macro-kinetics  with  a  dual-site  mechanism  for  the  Pt-Ru  porous 
anode  in  DMFC. 


2.3.  Generalizing  and  decoupling  the  model  equations 


Generalizing  the  model  equations  [7-11],  by  defining  dimen¬ 
sionless  variables  CM  =  cM/c^,  <2>  =  /3E,  X=x//,  the  dimensionless 
theoretical  expressions  and  dimensionless  parameters  can  be 
obtained  as  follows.  Moreover,  to  improve  the  model  and  the  pre¬ 
cision  of  the  solution,  the  boundary  condition  for  overpotential  at 
X=  1  is  normalized  by  defining  a  relative  dimensionless  overpoten¬ 
tial  ^  =  <PI/3E0  =  @I<P° ,  we  have 

j2  r 

^  =s-0oh,ru0co  exp  [PE0*), 
d2ds 

-jjjp  =  ^  '  ^oh.ru^co  exp (PE0^)  (18) 

with  boundary  conditions: 


X  =  1,  CM  =  1,  ^  =  1 


The  dimensionless  modulus  in  Eq.  (18),  s  =  ajol2/nFDec^  char¬ 
acterises  the  resistances  of  mass  transport  in  the  porous  anode. 


The  dimensionless  modulus  /z  =  v2/0E°  =  v2/<2>°,  where  v2  = 
ajol2  j3/ice,  characterises  the  relative  resistance  of  charge  transport 
when  applying  different  overpotentials  at  the  boundary  X=  1.  Thus, 
a  set  of  the  generalizing  model  Eq.  (21)  can  be  used  to  theoreti¬ 
cally  analyze  the  effects  of  the  coupled  activation,  concentration 
and  ohmic  polarisation  on  overall  polarisation  in  the  porous  anode 
of  the  DMFC. 

Decoupling  Eq.  (21 ),  we  have 


d2CM  _  s  d2xfs 
dX2  =  /z  Hx2 

With  the  boundary  conditions  we  obtain  from  Eq.  (19) 

(cm-i)=£op-i) 

/Z 

Substituting  Eq.  (20)  into  Eq.  (18),  we  have 


d2CM 

dX2 


=  S  •  0OH,Ru#CO  exp 


1  +  ^(CM  - 1 


} 


(19) 


(20) 


(21) 


X  =  0, 
X  =  1, 


dCM  _  0 
dX 

CM  =  1 


where  0Oh,ru  or  $co  is  a  function  of  CM  only. 

Solution  of  Eq.  (21)  thus  gives  the  concentration  distribution  in 
the  anode  and  from  Eq.  (20)  the  potential  distribution  also  can  be 
obtained. 


2.4.  Polarisation  curves  and  effectiveness  factors  of  the  porous 
anode 


According  to  Ohm’s  law,  the  local  current  density  at  a  lateral 
distance  x  is 

m  dcM  dE 

hoc  -  -nFDe-^-  -  -Kefa  (22) 

The  dimensionless  current  density  is  thus  defined  as 


I  nF 1  dCu 

Jk  =  “"fs¥ 


v2  dX 


2_  (Mf_ 

/z  dX 


(23) 


From  expressions  of  s  and  /z  (or  v2),  we  obtain  the  relation 
between  j  and  J,  i.e. 


l-ti-Xn,  1241 

Hence  Eqs.  (22)  and  (23)  describe  the  distributions  of  current 
density  and  the  total  current  density,  i.e.  the  apparent  (macro-) 
current  density  of  a  porous  electrode  with  a  thickness  l ,  will  be 


h  =  dioUj  =  -Ke 


—fiT  \~dX  ) x=i 


KeE°  fd^\ 

~T\dx)x= i 


(25) 


where  the  relation  of  E°-ji  describes  the  apparent  (macro-)  polar¬ 
isation  curve  of  the  porous  anode  and  the  relation  of  the 

apparent  dimensionless  polarisation  curve. 

The  corresponding  apparent  (macro-)  reaction  rate  can  be  writ¬ 
ten  as 


R=4  =  ^jT=_De 

nF  nFJ 


DeCM  (  dCM  \ 
l  V  dX  A=1 


(26) 


Introducing  the  effectiveness  factor  of  the  porous  anode  to  scale 
the  impact  of  parameters  such  as  a,  l,  De  and  xe, 


5  = 


apparent  rate  or  current  density 
intrinsic  rate  or  current  density 


h 

alKc°M,E°) 


(27) 


Table  2 

Physical  parameters  used  for  the  modeling. 
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Electrode  parameters 

References 

Catalyst  layer  thickness,  l  (m) 

3.8  x  10-6 

[4] 

5.0  x  10“6 

[32] 

1.0  x  10-5 

[19,29] 

2.3  x  10-5 

[18] 

2.5  x  10-5 

[30] 

5.0  x  10-5 

[13] 

Specific  area  of  anode,  a  (m-1  ) 

118,317 

[29] 

1.0  x  105 

[17] 

Porosity  of  anode,  e 

0.3 

[31] 

Diffusion  coefficients,  D0  (m2  s-1 ) 

2.8  x  10-9  exp[2436((l/353)  -  (1/T))] 

[20,31] 

1.6  x  10-9  (70°C) 

[18] 

6.03  x  10-9  (60  °C) 

[14] 

Effective  diffusion  coefficients,  De  (m2  s-1 ) 

De=Do£15 

[6,18] 

Effective  conductivity,  tce  (C2-1  m-1) 

3.4 

[19,29] 

Substituting  Eq.  (4)  into  Eq.  (27),  we  have 

Q(/o^oH,Ru^CO  exP(/^£°)  ^OH.Ru^CO  exP(^°) 

Hence  the  relationship  curves  of  jT-§  and  JT-§  can  be  calculated 
for  the  Pt-Ru  porous  anode. 

2.5.  Numerical  solution 

The  model  Eq.  (21)  was  solved  by  Newman’s  BAND  (j)  [6]  based 
on  a  finite  difference  method.  In  each  step,  AX/,  a  subroutine  for 
solving  a  set  of  nonlinear  algebraic  equations  was  inserted  to  cal¬ 
culate  the  values  of  #oh,ru(*')  and  6Co(i)  which  would  be  substituted 
into  the  source  term  [11].  Therefore  all  values  of  tf'(X)  and  CM(X) 
at  nodes  from  1  to  N+l  can  be  obtained;  namely  the  numerical 
solution  of  the  model  equation  [7,11  ].  Once  the  model  is  solved,  the 
polarisation  curve  as  well  as  effectiveness  factor  can  be  obtained 
from  the  grads  of  CM(1)  or  ^(1)  according  to  Eqs.  (25)  and  (28) 
[7,9]. 

3.  Analysis  of  the  macro-kinetics 

In  this  paper,  the  physical  parameters  used  are  taken  from 
the  literature  [4,6,13,14,17-20,29-32]  and  are  listed  in  Table  2. 
Furthermore  one  set  of  kinetic  parameters,  corresponding  to  the 
elementary  reaction  steps,  is  regressed  and  listed  in  Table  3,  which 


Table  3 

Kinetic  parameters  used  for  the  model. 


Kinetic  parameters 

30°C 

60  °C 

90°C 

fci  (ms-1) 

8.7  x  10-7 

4.2  x  10-6 

1.0  x  10-5 

k\  (molm-2  s-1) 

4.0  x  10-4 

1.5  x  10-3 

2.6  x  10-3 

/<2  (molm-2  s-1) 

3.5  x  10-9 

9.5  x  10-8 

8.0  x  10-7 

k 3,1  (molm-2  s-1) 

- 

3.2  x  10-9 

- 

k'3  x  (molm-2  s-1) 

- 

1.9  x  10-9 

- 

/<3,2  (mol  m-2  s-1) 

4.0  x  10-5 

5.0  x  10-5 

6.0  x  10-5 

k'3  2  (molm-2  s-1) 

3.0  x  10-5 

1.8  x  10-5 

1.4  x  10-5 

k 4  (molm-2  s-1) 

5.3  x  10-2 

5.9  x  10-2 

6.2  x  10-2 

0(2 

0.79 

ft  (V-1) 

0.5 

04  (V-1) 

0.5 

Fig.  2.  Variations  of  coverage  ratios  with  overpotential  at  different  temperatures. 


fit  with  the  experimental  polarisation  data  for  bulk  methanol  con¬ 
centrations  0.5  M,  1.0 M  and  2.0  M  (mol  dm-3)  and  temperatures  of 
30  °C,  60  °C  and  90  °C,  reported  by  Scott  in  Ref.  [15]. 

3.1  Coverage  ratios  of  species  on  the  electrode  surface 

3.1.1.  Variations  of  coverage  ratios  with  overpotential 

Generally  the  superficial  rate  or  current  density  is  proportional 
to  both  coverage  ratios  0Oh,ru  and  6C o,  which  are  mutually  coupled 
with  the  coverage  ratio  0m  of  CH3OHads  adsorbed  at  Pt  sites.  Also  all 
of  the  coverage  ratios  are  functions  of  concentration,  temperature 
and  overpotential. 

Fig.  2  shows  the  variations  of  surface  coverage  ratios  with  over¬ 
potentials  E,  calculated  by  using  the  kinetics  parameters  at  30  °C, 
60  °C  and  90  °C  with  1.0  M  of  methanol  for  Case  A.  The  values  for 
6m  and  0Oh,ru  are  shown  on  the  left  y-axis  while  values  for  6C o  are 
shown  on  the  right  y-axis  with  a  scale  of  10-5. 

With  increasing  values  of  E,  $oh,ru  rapidly  increases  to  near  1, 
undergoes  a  slight  decrease  and,  approaches  1  again:  this  means  the 
electrochemical  formation  of  OH  at  Ru  site  (step  3)  may  not  become 
a  rate  limiting  step.  6M  initially  does  not  change  with  an  increase  in 
potential  (maintains  a  value  approximately  that  of  the  adsorption 
equilibrium)  but  then  undergoes  a  sharp  decrease,  finally  maintain¬ 
ing  a  small  near  constant  value.  This  means  that  the  adsorption  rate 
of  methanol  (step  1 )  may  be  the  rate  limiting  step  with  high  overpo¬ 
tentials.  However,  values  for  6C o  are  always  small;  this  means  that 
the  rate  of  surface  reaction  (last  step  4)  is  much  quicker  than  the  rate 
of  step  2.  With  increasing  potential,  6C o  first  shows  a  sharp  increase 
and  then  undergoes  a  sharp  decrease,  ultimately  maintaining  very 
small  values;  this  means  that  electrochemical  dehydrogenation  of 
CH3OHads  at  Pt  sites  (step  2)  may  be  the  rate  limiting  step  at  low 
overpotentials.  Hence  there  is  no  rate  limiting  step  in  the  transition 
region  with  intermediate  values  of  E. 

In  the  model  the  influence  of  temperature  on  the  coverage  ratios 
is  important.  The  activation  energies,  determined  from  Arrhenius 
plots  of  the  rate  constants  (Table  2)  with  inverse  temperature  (not 
shown),  of  methanol  adsorption  of  step  1  and  of  CH3OHads  dehydro¬ 
genation  of  step  2  were  approximately  40  kj  mol-1  and  80  kj  mol-1 
respectively.  Looking  back  to  Fig.  2,  we  can  see  that  the  greatest 
influence  of  temperature  on  the  coverage  ratios  will  take  place  in 
the  case  without  a  limiting  step  and  that  0C o  is  more  sensitive  to 
temperature  than  the  other  two  coverage  ratios;  especially  in  the 
transition  region.  This  indicates  that  the  apparent  activation  energy 
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Fig.  3.  Variations  of  coverage  ratios  with  overpotential  at  different  methanol  con¬ 
centrations. 

in  the  transition  region  should  be  higher  than  that  of  any  step  of  the 
series  of  elementary  reactions;  its  maximum  is  equal  to  the  sum  of 
activation  energies  of  all  consecutive  reactions,  according  to  the 
kinetics  of  multistep  reactions  [22]. 

Fig.  3  shows  the  variation  in  surface  coverages  with  potential  at 
different  concentrations  for  Case  A  at  60  °C.  The  values  of  0m  appear 
to  be  more  sensitive,  and  the  values  of  0Oh,ru  and  0C o,  less  sensitive 
to  concentration  than  to  temperature.  Similar  to  the  trends  in  the 
data  shown  in  Fig.  2,  the  rate  limiting  step  will  be  step  2  at  low  E, 
then  changes  to  step  1  at  high  E  with  increasing  overpotential. 

The  current  density  of  the  limiting  step  1  can  be  described 
approximately  by  a  kinetic  expression  j  =  nFk^c^  (j  is  proportional 
to  the  bulk  c^)  as  the  values  of  0m  and  0C o  approach  zero  at  high  E , 
according  to  the  rate  expression  for  r\  in  Table  1  for  Case  A. 

Theoretically,  the  last  step,  4,  of  surface  reaction  producing  C02 
might  be  a  rate  limiting  step  under  certain  conditions  such  as  high 
temperatures  and  concentrations,  but  not  with  the  kinetic  values 
adopted  in  this  paper. 

Overall,  the  rate  limiting  step  may  change  from  step  2,  indepen¬ 
dent  of  concentration,  and  step  1,  independent  of  overpotential, 
which  will  not  only  lead  to  a  change  in  apparent  activation  energy 
but  also  to  a  change  of  reaction  order  between  0  and  1  with  respect 
to  methanol  concentration.  Furthermore  it  should  be  noted  that 
Oqq  cannot  be  neglected  since  it  is  proportional  to  the  total  reaction 
rate,  although  its  values  are  quite  small.  The  sharp  change  of  0C o 
as  well  as  0M  with  overpotential  is  the  reason  that  the  kinetics  of 
methanol  oxidation  on  Pt-Ru  cannot  be  simply  described  by  B-V 
equations. 

3.1.2.  Distributions  of  coverage  ratios  in  the  porous  anode 

Spatial  distributions  in  coverage  ratios  will  exist  in  the  porous 
anode,  as  the  coverage  ratios  are  functions  of  concentration,  tem¬ 
perature  and  overpotential.  Flence  these  distributions  will  also 
influence  the  local  current  density  (or  local  reaction  rate)  dis¬ 
tribution  in  the  porous  anode.  Since  the  product  ^oh.ru^co  is 
directly  proportional  to  surface  current  density,  the  distributions 
of  ^oh.ru^co  will  have  considerable  influence  on  the  local  current 
density. 

The  variations  in  the  relative  local  current  density  jloc(x)/jloc(0) 
distributions  with  the  dimensionless  overpotential  0°  are  shown 
in  Fig.  4.  This  example  was  calculated  by  solving  the  macro-kinetic 


Fig.  4.  Variation  of  distribution  of  relative  j\oc  with  0°. 


model  equations  at  pore  depth  x  =  0(xm,  lOpm,  25  pan,  50  pan 
respectively  for  a  porous  electrode  of  thickness  l  =  50  pm.  From  the 
relative  values  of  jloc  to  the  superficial  jloc(0),  it  can  be  seen  that 
jloc(x)  is  evidently  impaired  by  coupled  mass  and  charge  transport 
for  longer  pores,  at  higher  E. 

A  simpler  way  to  solve  the  model  equations  for  macro-kinetics 
is  to  take  the  product  (^oh,ru^co)/=o  as  independent  of  x  and  not  a 
function  of  x;  however  this  will  bring  considerable  errors.  By  check¬ 
ing  computations,  the  relative  errors  of  local  current  densities  j\oc 
(^oh,ru^co)/-x/Jioc  (#oh,ru#co)i-o  are.  for  example,  from  20%  to  50%, 
when,  x  >  10  |xm  and  <J>°  >  6.0.  Therefore  it  is  important  to  involve 
the  coverage  ratio  distributions  in  the  analysis  of  macro-kinetics  of 
the  porous  anode  at  most  overpotentials. 

3.2.  Macro-polarisation  curve  and  effectiveness 

3.2.1.  Macro -polarisation  curve 

The  polarisation  curve  of  a  porous  electrode,  called  the  macro¬ 
polarisation  curve,  describes  the  variation  in  the  applied  surface 
overpotential  E°  and  integral,  total  current  density  jj. 

Fig.  5a  shows  examples  of  the  dimensional  macro-polarisation 
curves  corresponding  to  physical  parameters  arranged  into  five 
groups  by  3  levels  of  electrode  thickness,  /,  2  levels  of  specific  area, 
a  and  1  level  of  effective  diffusion  coefficient,  De.  The  parameters 
used  for  the  curve  numbered  4  (Fig.  5a)  are  very  close  to  those 
used  in  Scott’s  papers  [15,21].  All  the  macro-polarisation  curves 
possess  the  same  trend;  an  approximate  exponential  rise  in  current 
density  with  increasing  overpotential  and  eventually  approaching 
a  maximum  current  density.  This  maximum  current  density  is  an 
“adsorption  limiting  current  density”,  ja\  rather  than  the  familiar 
“limiting  current  density”  caused  by  external  mass  transport. 

The  data  show  that  thicker  electrodes  and  higher  specific  elec¬ 
trode  areas  result  in  higher  current  densities  at  a  given  potential. 
The  current  densities  are  determined  by  the  very  low  dehydrogena¬ 
tion  rate  of  step  2  at  low  overpotentials  as  well  as  by  the  adsorption 
rate  of  step  1  at  high  overpotentials. 

Fig.  5b  shows  the  dimensionless  macro-polarisation  curves J-0° 
corresponding  to  the  same  physical  parameters  of  Fig.  5a,  and  thus 
the  dimensionless  values  of  v2  and  s.J  is  proportional  to  the  mean 
value  of  current  density  per  unit  of  inner  surface  area  over  the 
porous  anode  jT/a/,  and  therefore  the  J-0°  curve  shows  the  effects 
of  mass  and  charge  transfer  on  the  porous  anode.  The  curve  with 
symbol  A  represents  the  dimensionless  intrinsic  polarisation  curve 
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Fig.  5.  (a)  Macro-polarisation  curves  of  jT-£°  for  different  a  and  /.  (b)  Dimensionless 
macro-polarisation  curves  J-<P°  for  different  v2  and  s. 


independent  of  physical  parameters,  and  can  be  calculated  by  the 
macro-model  equation,  but  with  1^0,  say  /  =  10— 8  m,  i.e.  both  v2 
and  s  ->  0  in  the  numerical  procedure. 

The  decreasing  values  of  v2  and  s  corresponding  to  curves  of 
1,  2,  3,  4  and  5  result  in  higher  current  densities  and  indicate  the 
decreasing  influence  of  transport  processes  in  the  porous  electrode. 
High  values  of  v2  and  s  (approximately  40-50)  for  curve  1  have  a 
large  effect  of  transport  while  low  values  (<1.0)  for  curves  4  or  5 
near  the  intrinsic  curve  (denoted  by  symbol  A),  had  a  slight  effect. 
The  physical  parameters  used  for  curve  4  are  close  to  those  in  Scott’s 
papers  [15,21  ],  which  means  that  Scott’s  kinetic  experiments  were 
carried  out  under  approximate  intrinsic  conditions. 

3.2.2.  Effectiveness  of  the  porous  anode 

The  effects  of  mass  and  charge  transport  can  be  characterised 
by  the  effectiveness  §  of  the  porous  anode.  The  %-J  curves  calcu¬ 
lated  corresponding  to  the  physical  parameters  in  Fig.  5,  are  shown 
in  Fig.  6a.  The  effectiveness  is  higher  for  thinner  electrodes  and  for 
higher  specific  areas  due  to  the  decreasing  influence  of  mass  trans¬ 


Fig.  6.  (a)  Variation  of  efficiency  factor  f  with  dimensionless  current  density  J  for 
different  /  and  a.  (b)  Variation  of  effectiveness  £  with  dimensionless  current  density 
J  for  different  l  and  De. 


port  processes.  With  an  increase  in  macro-current  density,  every 
point  of  the  porous  anode  can  undergo  a  change  of  limiting  step 
from  step  2  to  step  1,  depending  on  the  local  concentration  and 
overpotential  at  the  point.  Therefore,  the  end  point  of  each  curve 
represents  the  dimensionless  adsorption  limiting  current  density, 
Jai,  where  only  the  effect  of  mass  transport  exists  since  the  limiting 
step  1  is  independent  of  overpotential.  The  data  with  symbol  □  over¬ 
lapping  curve  4  (Fig.  6a)  are  calculated  from  the  physical  parameters 
used  in  Scott’s  papers  [15,21  ] :  all  have  values  of  §  greater  than  0.96, 
approaching  a  value  of  1.0. 

Theoretically  the  effectiveness  is  a  function  of  both  the  effec¬ 
tive  diffusion  coefficient  De  and  conductivity /ce  characterizing  both 
mass  and  charge  transport  in  the  porous  anode.  Fig.  6b  shows  the 
variations  of  §  with  De  which  essentially  considers  the  effect  of 
porosity  £  of  the  anode,  through  for  example  the  Bruggeman  equa¬ 
tion  (i.e.  De  =  D0£1-5)  [6,18].  The  physical  parameters  are  arranged 
into  nine  groups  by  three  levels  of  electrode  thickness,  l,  three 
levels  of  effective  diffusion  coefficient  De  and  one  level  of  spe¬ 
cific  area  a.  The  three  groups  of  curves,  for  the  same  value  of  De 
show  the  effects  of  e  and  l  on  §.  For  a  thinner  porous  electrode 
as  well  as  a  larger  De,  the  influence  of  pore  diffusion  becomes 
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less  significant,  i.e.  higher  effectiveness  values  are  achieved.  Gener¬ 
ally,  effectiveness  decreases  with  increasing  macro-current  density. 
However,  for  a  thicker  porous  electrode  as  well  as  a  larger  De,  effec¬ 
tiveness  goes  through  a  minimum  with  current  density,  which  is 
caused  by  the  effects  of  both  mass  and  charge  transport.  With 
unceasingly  increasing  current  density,  the  effect  of  charge  trans¬ 
port  resistance  gradually  disappears  due  to  the  change  in  the 
limiting  reaction  step  (methanol  adsorption).  The  effect  of  bet¬ 
ter  mass  transport,  due  to  the  larger  De,  means  that  more  of  the 
electrode  would  be  operating  at  or  close  to  a  limiting  current 
and  hence  effectiveness  increases  again  as  J  increases,  until  Jal  is 
reached. 

4.  Verification  of  the  macro-kinetic  model 

4  A.  Type  of  reaction  mechanism 

For  model  verification,  the  predicted  macro-polarisation  curves 
were  compared  to  the  experimental  polarisation  data  reported  by 
Scott  or  Nordlund’s  papers  for  methanol  oxidation  on  the  Pt-Ru 
porous  anodes  [15,16,18]. 

Fig.  7  shows  two  predicted  macro-polarisation  curves  (solid 
lines  1  and  2),  and  two  intrinsic  polarisation  curves  (dashed  lines  V 
and  2'),  for  the  two  kinetic  models  Case  A  (noted  1  and  1')  and  Case 
B  (noted  2  and  2')  respectively.  The  intrinsic  polarisation  curves 
were  calculated  by  Eq.  (3)  j  =  nFl-r^c^,  E0)]  and  are  independent 
of  any  physical  parameters,  which  also  consistent  with  those  calcu¬ 
lated  by  the  macro-model  equation  but  as  the  thickness  approached 
the  limit,  l  ->  0,  in  the  numerical  procedure.  The  experimental  data, 
which  are  close  to  the  model  data  for  the  intrinsic  situation,  from 
Scott’s  paper  [15]  are  denoted  by  the  symbol  Q.  The  four  curves 
only  separate  at  high  current  densities  due  to  the  effects  of  trans¬ 
ports  or  different  mechanisms.  The  influence  of  the  mechanism 
can  be  appreciated  (see  Fig.  8)  from  the  differences  in  the  varia¬ 
tion  surface  coverages  for  the  two  mechanisms.  Fig.  8  shows  for 
Case  B,  the  variation  in  coverage  ratios  with  overpotential,  using 
the  kinetics  parameters  at  60° C  at  1.0  M  methanol.  The  main  differ¬ 
ence  between  the  curves  of  Case  B  and  Case  A  (Figs.  8  and  2)  is  that 
the  coverage  ratio  0oH,pt  sharply  increases  at  high  overpotentials: 
this  makes  the  polarisation  curves  turn  back  (decrease  in  current 
density)  at  high  potentials  due  to  the  loss  of  Pt  vacancy  sites  for 
methanol  adsorption  [15,18,23]. 


Fig.  7.  Polarisation  curves  of  jj-E°  for  both  Case  A  and  Case  B  mechanisms. 


Fig.  8.  Variations  of  coverage  ratios  (0M,  0oh,ru.  0c o  and  0OH,pt)  with  overpotential  of 
mechanism  Case  B. 


4.2.  Effect  of  temperature  and  bulk  concentration 

Fig.  9a  shows  the  predicted  intrinsic  (dashed  lines)  and  macro¬ 
polarisation  curves  (solid  lines)  using  at  three  temperatures:  30  °C, 
60  °C  and  90  °C.  The  experimental  data  from  Scott  and  co-workers 
[  1 5  ]  are  also  shown  for  different  temperatures.  As  expected  increas¬ 
ing  temperature  increases  current  density.  The  data  show  again  that 
Scott’s  kinetic  experiments  were  carried  out  under  approximate 
intrinsic  conditions.  The  predicted  values  agree  well  with  experi¬ 
mental  data  at  lower  currents  at  30  °C  and  60  °C,  although  there  is 
a  slight  departure  from  the  data  at  higher  currents,  which  shows  a 
linear  trend,  at  90  °C. 

Fig.  9b  shows  the  predicted  intrinsic  (dished  lines)  and  macro¬ 
polarisation  curves  (solid  lines)  at  60  °C  for  three  concentrations  of 
methanol  0.5  M,  1.0  M  and  2.0  M.  The  experimental  data  from  Scott 
and  co-workers  [15]  are  indicated  for  different  concentrations.  The 
predicted  values  agree  well  with  experimental  data  at  lower  cur¬ 
rent  densities  for  0.5  M  and  1.0  M  concentrations;  however  there  is 
a  considerable  departure  from  the  data  at  higher  current  densities 
at  2.0  M.  The  experimental  current  densities  for  a  concentration  of 
2.0  M  are  lower  than  those  predicted,  and  close  to  the  data  with 
cjjj  =  1 .0  M.  A  possible  explanation  for  this  may  be  that  the  effect  of 
mass  transport  on  the  exterior  surface  of  the  porous  anode  becomes 
important  at  high  currents  as  C02  bubbles  may  block  methanol 
transfer  from  the  bulk  into  the  pores  of  the  anode. 

4.3.  Effect  of  anode  activity  and  thickness 

Fig.  10a  and  b  shows  the  predicted  intrinsic  (dished  lines)  and 
macro-polarisation  curves  (solid  lines)  for  a  porous  anode  with 
a  higher  activity  and/or  larger  thickness.  There  are  considerable 
differences  between  the  intrinsic  and  macro-polarisation  curves 
due  to  the  influence  of  transport  processes  in  the  anode.  Generally 
transport  limitation  decreases  the  current  density  from  the  anode 
at  the  same  overpotential. 

In  Fig.  10a  and  b  the  experimental  polarisation  data  reported 
by  Nordlund  and  Lindbergh  are  shown  [14,18].  Porous  anodes  pre¬ 
pared  by  different  research  groups  may  have  different  activities, 
microstructures,  etc.  The  thickness  (1  =  23  pan)  of  the  porous  anode 
used  to  obtain  the  data  [18]  in  Fig.  10b  is  approximately  twice  that 
(/  =  14  fun)  for  the  data  in  Fig.  10a  [14].  Moreover  the  experimental 
current  densities  in  Fig.  10a  and  b  are  greater  than  those  of  Scott 
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Fig.  9.  Polarisation  curves  of  jT-E°  for  (a)  different  temperatures  and  (b)  bulk  con¬ 
centrations  of  methanol. 

and  co-workers  [15]  in  Fig.  9a  and  b.  Hence  we  modified  the  values 
of  kinetic  parameters  which  are  listed  in  Table  4  for  the  calculation 
of  intrinsic  and  macro-polarisation  curves  at  70  °C. 

In  Fig.  10a,  the  data  with  cj^  of  0.5  M  (denoted  by  symbol  □), 
is  consistent  with  the  predicted  macro-polarisation  curve;  how¬ 
ever  the  data  with  of  1.0  and  2.0  M  (denoted  by  symbols  O 
and  A)  evidently  depart  from  the  corresponding  predicted  macro¬ 
polarisation  curves  at  higher  currents,  i.e.  the  linear  trend  of  the 
data.  In  Fig.  10b,  for  the  porous  electrode  with  a  higher  activity 
as  well  as  a  larger  thickness,  the  data  with  of  1.0  M  and  2.0  M 
are  still  consistent  with  the  predicted  macro-polarisation  curves  at 

Table  4 

Kinetic  parameters  based  on  Ref.  [15]  at  70  °C. 


Fig.  10.  Polarisation  curves  of/T-E°  for  the  anode  with  a  (a)  higher  activity  and  (b) 
higher  activity  and  a  larger  thickness. 

lower  currents,  but  depart  considerably  from  the  predicted  macro¬ 
polarisation  curves  at  higher  currents.  The  data  for  2.0  M  methanol 
becomes  much  closer  to  that  for  1.0  M  and  actually  lie  on  the  pre¬ 
dicted  macro-polarisation  curve  for  1.0  M. 

As  mentioned  above,  the  macro-kinetic  model  developed  here 
does  not  include  the  influence  of  two-phase  flow  caused  by  C02 
bubbles,  which  impede  methanol  transport.  Thus  is  likely  that  the 
departure  of  experimental  data  from  the  predicted  data,  tending 
to  be  linear  at  higher  currents,  may  be  that  the  electrode  polari¬ 
sation  here  is  dominated  by  concentration  polarisation  caused  by 
mass  transfer  limitations  of  two-phase  flow  rather  than  the  kinetic 
polarisation  at  higher  currents  as  more  C02  is  produced. 

5.  Conclusion 

A  macro-kinetic  model  of  a  porous  Pt-Ru  anode  for  methanol 
oxidation  has  been  developed.  The  performance  of  the  porous 
Pt-Ru  anode  depends  on  the  mechanism  and  parameters  of  the 
intrinsic  kinetics  of  methanol  oxidation  as  well  as  parameters  such 
as  thickness  of  the  electrode,  surface  area,  effective  diffusion,  charge 
transfer  coefficients,  etc. 
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Intrinsic  kinetics  of  the  oxidation  are  only  related  to  the  flux  of 
charge  on  the  surface  of  porous  electrode,  i.e.  the  current  density 
being  independent  of  any  physical  parameter,  but  macro-kinetics 
are  related  to  the  divergence  of  current  density  in  the  volume  of  the 
porous  electrode,  the  latter  requires  a  set  of  second  order  differen¬ 
tial  equations  to  describe  the  coupled  physical  transport  process 
and  chemical  reaction  processes  in  a  three  dimension  space. 

The  mathematical  model  of  macro-kinetics,  developed,  has  two 
traits:  firstly  it  contains  the  dual-site  mechanism  of  methanol 
oxidation  on  Pt-Ru  anode;  secondly  the  model  equations  are  gen¬ 
eralized. 

By  using  the  finite  difference  method  with  a  subroutine  for  solv¬ 
ing  a  set  of  nonlinear  algebraic  equations  in  each  step,  the  solutions 
of  the  model,  the  distributions  of  concentration  and  overpotential, 
were  obtained  as  well  as  the  surface  coverage  ratio  distribu¬ 
tion,  effectiveness  of  the  porous  electrode  and  macro-polarisation 
curves. 

Comparing  the  predicted  polarisation  curves  with  experimental 
polarisation  data,  the  theoretical  analysis  for  the  performance  of  the 
porous  anode  can  be  carried  out  to  understand  macro-kinetics  as 
well  as  to  find  out  ways  to  improve  the  performance  of  the  electrode. 

From  the  theoretical  analysis,  the  apparent  current  densities  are 
very  small  and  almost  approach  zero  at  low  overpotentials  for  all 
cases  since  the  dehydrogenation  rates  of  CH3OHads  adsorbed  on 
the  Pt  sites  are  too  small.  Therefore  improving  the  electrocatalyst 
activity  is  an  important  issue  [24,25]  for  reducing  the  considerable 
electrode  polarisation  near  “open-circuit”.  Comparing  the  kinetic 
parameters  in  Table  3  with  those  in  Table  2,  it  seems  that  the  capac¬ 
ity  of  the  Pt-Ru  anode  for  methanol  oxidation  is  actually  limited  by 
the  rate  of  methanol  adsorption  on  the  Pt  sites;  hence  to  accelerate 
the  adsorption  on  the  electrocatalyst  is  another  important  issue. 

The  predicted  polarisation  curves  calculated  by  the  model  of 
macro-kinetics  without  taking  into  account  two-phase  flow,  are 
consistent  with  experimental  polarisation  data  at  lower  current 
densities.  This  feature  implies  that  the  porous  electrode,  with  a 
high  specific  area  for  higher  apparent  current  density  at  lower  local 
current  density,  will  be  instrumental  in  reducing  the  influence  of 
two-phase  flow.  The  departure  of  experimental  data  from  the  pre¬ 
dicted  polarisation  curves  at  high  concentration  at  high  apparent 
current  densities  is  likely  due  to  the  two-phase  flow  mode.  Opti¬ 
mising  the  electrode  structure  in  relation  to  the  two-phase  flow  is 
also  an  important  issue  for  further  investigation  [26-29]. 

This  work  also  can  be  referenced  for  theoretical  analyses  of 
porous  or  packed  bed  anode  in  other  Direct  Liquid-Feed  Fuel  Cell, 
using  other  alcohols,  borohydride,  etc. 


Acknowledgement 

The  authors  acknowledge  the  support  of  research  project  NSFC 
(20776091 )  provided  by  China  National  Natural  Science  foundation. 

References 

[1  ]  K.  Scott,  A.  Shukla,  in:  R.E.  White  (Ed.),  Modern  Aspects  of  Electrochemistry,  No. 
40,  Springer,  New  York,  2007,  pp.  127-228. 

[2]  B.  Garcia,  J.  Weidner,  in:  R.E.  White  (Ed.),  Modern  Aspects  of  Electrochemistry, 
No.  40,  Springer,  New  York,  2007,  pp.  229-284. 

[3]  S.  Gottesfeld,  J.  Power  Sources  171  (2007)  37-45. 

[4]  J.  Meyers,  J.  Newman,  J.  Electrochem.  Soc.  149  (6)  (2002)  A718-A728. 

[5]  J.  Newman,  C.  Tobias,  J.  Electrochem.  Soc.  109  (12)  (1962)  1183-1191. 

[6]  J.  Newman,  Electrochemical  Systems,  third  ed.,  Prentice-Hall,  Newjersey,  2004. 

[7]  K.  Scott,  Y.-P.  Sun,  in:  C.  Vayenas,  M.  Gamboa-Adelco  (Eds.),  Modern  Aspects  of 
Electrochemistry,  No.  41,  Springer,  New  York,  2007,  pp.  221-304. 

[8]  Y.-P.  Sun,  K.  Scott,  Fuel  Cell  4  (1/2)  (2004)  30-38. 

[9]  Y.-P.  Sun,  K.  Scott,  Chem.  Eng.  J.  1024  (1)  (2004)  83-91. 

[10]  Y.-P.  Sun,  J.  Chem.  Ind.  Eng.  58  (9)  (2007)  2161-2168  (in  Chinese). 

[11]  Y.-P.  Sun,  L.  Xing,  CIESC  60  (1)  (2009)  55-68  (in  Chinese). 

[12]  H.  Gasteiger,  N.  Markovic,  P.  Ross  Jr.,  J.  Cairns,  J.  Phys.  Chem.  97  (1993) 
12020-12029. 

[13]  P.  Kauranen,  E.  Skou,  J.  Munk,  J.  Electroanal.  Chem.  404  (1996)  1-13. 

[14]  J.  Nordlund,  G.  Lindbergh,  J.  Electrochem.  Soc.  151  (9)  (2004)  A1357-A1362. 

[15]  M.  Shivhare,  R.  Allen,  K.  Scott,  A.  Morris,  E.  Martin,  J.  Electroanal.  Chem.  595 
(2006)  145-151. 

[16]  M.  Shivhare,  C.  Jackson,  K.  Scott,  E.  Martin,  J.  Power  Sources  173  (2007)  240- 
248. 

[17]  B.  Garcia,  V.  Sethuraman,  J.  Weidner,  R.E.  White,  R.  Dougal,  J.  Fuel  Cell  Sci. 
Technol.  1  (2004)  43-48. 

[18]  J.  Nordlund,  G.  Lindbergh,  J.  Electrochem.  Soc.  149  (9)  (2002)  A1107-A1113. 

[19]  K.  Scott,  P.  Argyropoulos,  J.  Power  Sources  137  (2004)  228-238. 

[20]  K.  Scott,  W.  Taama,  S.  Kramer,  P.  Argyropoulos,  K.  Sundmacher,  J.  Electrochim. 
Acta  45  (1999)  945-957. 

[21]  I<.  Scott,  P.  Argyropoulos,  K.  Sundmacher,  J.  Electroanal.  Chem.  477  (1999) 
97-110. 

[22]  F.G.  Helfferich,  Kinetics  of  Multistep  Reactions,  Elsevier,  Amsterdam,  2004. 

[23]  P.  Kauranen,  E.  Skou,  J.  Munk,  J.  Electroanal.  Chem.  404  (1996). 

[24]  N.  Tian,  Z.-Y.  Zhou,  S.-G.  Sun,  Y.  Ding,  Z.-L.  Wang,  Science  316  (5825)  (2007) 
732-735. 

[25]  Q.  Wang,  G.Q,  Sun,  L.  Gao,  L.H.  Jiang,  G.X.  Wang,  S.L.  Wang,  S.H.  Yang,  Q.  Xin,  J. 
Power  Sources  177  (2008)  142-147. 

[26]  J.  Nordlund,  C.  Picard,  E.  Birgersson,  M.  Vynnycky,  G.  Lindbergh,  J.  Appl.  Elec¬ 
trochem.  34  (2004)  763-770. 

[27]  R.  Allen,  L.  Chan,  L.-X.  Yang,  K.  Scott,  S.  Roy,  J.  Power  Sources  143  (2005)  142- 
149. 

[28]  R.  Chetty,  K.  Scott,  J.  Electrochim.  Acta  52  (2007)  4073-4081. 

[29]  K.  Scott,  P.  Argyropoulos,  J.  Electroanal.  Chem.  567  (2004)  103-109. 

[30]  L.  Chan,  R.G.  Allen,  K.  Scott,  J.  Power  Sources  161  (2006)  11-18. 

[31  ]  K.  Scott,  W.  Taama,  J.  Cruickshank,  J.  Power  Sources  65  (1997)  159-171. 

[32]  Z.G.  Shao,  W.F.  Lin,  F.Y.  Zhu,  P.A.  Christensen,  M.Q,  Li,  H.M.  Zhang,  J.  Electrochem. 
Commun.  8  (2006)  5-8. 


